Mechanical Systems that are both Classical and Quantum 
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Abstract 

Quantum dynamics can be regarded as a generaliza- 
tion of classical finite-state dynamics. This is a fa- 
miliar viewpoint for workers in quantum computa- 
tion, which encompasses classical computation as a 
special case. Here this viewpoint is extended to me- 
chanics, where classical dynamics has traditionally 
been viewed as a macroscopic approximation of quan- 
tum behavior, not as a special case. 

When a classical dynamics is recast as a special 
case of quantum dynamics, the quantum description 
can be interpreted classically. For example, some- 
times extra information is added to the classical state 
in order to construct the quantum description. This 
extra information is then eliminated by representing 
it in a superposition as if it were unknown informa- 
tion about a classical statistical ensemble. This usage 
of superposition leads to the appearance of Fermions 
in the quantum description of classical lattice-gas dy- 
namics and turns continuous-space descriptions of 
finite-state systems into illustrations of classical sam- 
pling theory. 

A direct mapping of classical systems onto quan- 
tum systems also allows us to determine the mini- 
mum possible energy scale for a classical dynamics, 
based on a localized rate of state change. We use a 
partitioning description of dynamics to define local- 
ity, and discuss the ideal energy of two model systems. 



1 Introduction 

A digital photograph looks continuous but in fact, 
if you examine it closely enough, you discover that 
there is only a finite amount of detail. Similarly, a 



*MIT, Cambridge MA 02139; nhm@mit.edu 



digital movie looks like it is changing continuously in 
time, but in fact it is actually a discrete sequence of 
digital images. 

Something similar is true of nature. Although 
the world looks to our senses as if it has an infinite 
amount of resolution in both space and time, in fact 
a finite-sized physical system with a finite energy has 
only a finite amount of distinguishable detail and this 
detail changes at only a finite rate [?f] . 

1.1 A bit of history 

The finite character of the states of physical systems 
came as a great surprise when it became apparent 
at the start of the twentieth century. The revolution 
was started by Max Planck in 1900 when he found 
that he had to introduce a new constant into physics 
in order to understand the thermodynamics of elec- 
tromagnetic radiation in a cavity. The new constant 
fixed the statistical mechanical analysis, but it did so 
by making the count of distinct possible states finite. 

Planck's constant has a particularly simple inter- 
pretation in classical statistical mechanics. There, 
the number of microscopic ways a system can real- 
ize a macroscopic state is taken to be proportional to 
the volume in phase space (i.e., position/momentum 
space) of states consistent with the parameters that 
define the macroscopic state. In units where Planck's 
constant is one, the phase space volume becomes the 
actual finite count of distinct states. Planck's con- 
stant is the fundamental grain size in phase space, 
representing one state. 

Although Planck's constant was initially intro- 
duced to fix statistical mechanics, a new dynami- 
cal theory also grew from this beginning. The re- 
sulting quantum mechanics established new rules for 
describing the microscopic dynamics of physical sys- 
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terns. The QM rules were, of course, constructed so 
that the well established rules of classical mechanics 
would be recovered as a limiting case. The continuous 
equations of motion in CM correspond to the limit 
in which Planck's "fundamental grain size in phase 
space" goes to zero. In this limit, a finite-sized phys- 
ical system has an infinite number of distinguishable 
states and passes through an infinite subset of these 
states in a finite amount of time. Even though an 
actual finite-sized physical system has only a finite 
number of distinct states and a finite rate of state 
change, if these finite numbers are large enough (as 
they would normally be in a macroscopic system) the 
motion is well approximated by the continuous CM 
equations. 

1.2 Classical finite-state models 

The finite-state and finite-rate character of physical 
systems make classical finite state models of funda- 
mental interest in physics. This is because classical 
models are easier to understand and simulate than 
quantum models. 

In statistical mechanics, classical finite-state lat- 
tice systems have a long history as models of phase 
change phenomena in magnetic materials and mix- 
tures of fluids [32 [THl HH |5] . Such classical lattice 
gases can be considered special cases of general quan- 
tum lattice systems and described and analyzed us- 
ing the same quantum formalism |33j . By looking at 
these special cases, we gain insight into macroscopic 
classical concepts such as entropy and phase change 
by seeing them arise from underlying classical finite- 
state systems. 

We might hope to similarly gain insight by study- 
ing classical special cases of quantum rfynamicsj^ 
Any such classical special case must (like any other 
finite physical system) necessarily have only a finite 
number of distinct states. 

The most common classical finite-state models of 
physical dynamics are numerical approximations of 
continuum models. These are not good candidates 
to be special cases of microscopic physical dynamics. 



^Although the goal is very different, some work that seeks 
to describe all physical dynamics as classical also focuses on 
situations where QM acts classically [TSl [ill l36l l38l |42] . 



though, because in this case it is the theoretical con- 
tinuum model that has realistic physical character- 
istics, not its truncated and rounded-off finite-state 
realization [59] . The best candidates for exact special 
cases are lattice dynamics generalizations of classi- 
cal lattice gases. These are lattice models with exact 
finite-state conservation laws, that at large scales can 
exhibit not only realistic thermodynamic behavior 
but also realistic classical mechanical behavior such 
as wave propagation and fiuid dynamics [3U E5]j^ 
These models have been used to simulate physical 
dynamics [3S] but have not previously been studied 
as classical special cases of quantum dynamics. 

1.3 Finite-state physical dynamics 

Any invertible classical logic operation can be im- 
plemented as a QM unitary dynamics, and so can 
any composition of such operations. Since there are 
additional operations needed to implement an arbi- 
trary QM unitary dynamics [1], classical finite-state 
dynamics is logically a proper subset of quantum dy- 
namics. 

Before the advent of invertible models of computa- 
tion [Kig and QM models of computation [51IM1I5] it 
wasn't obvious that, from a logical standpoint, a clas- 
sical dynamics could be regarded as a special case of 
quantum dynamics. Historically, the viewpoint had 
been that classical dynamics is a purely macroscopic 
concept, relevant only in a limit in which quantum 
effects become negligible. In retrospect the confusion 
stems from not considering the possibility of invert- 
ible classical finite-state dynamics — an understand- 
able omission since QM was developed before com- 
puter ideas were prevalent, let alone the idea of in- 
vertible computation! 

We are glossing over a significant issue, though, 
when we picture classical finite-state dynamics as a 
composition of unitary operations. In a quantum 
computer [TJ [S] , the structures that switch from one 
unitary operation to the next are assumed to be 
macroscopic and classical. Thus it seems that our pic- 

^Discrcte classical models of general relativity [15[. devel- 
oped for use in quantizing gravity, might be regarded as a step 
towards describing gravity as a classical special case of finite- 
state physical dynamics. 
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ture is incomplete since there are unspecified things 
going on to implement the explicitly time-dependent 
sequence of unitary operations. It is helpful to note, 
though, that if we take the limit of an infinitely 
large array of finite-state elements going through a 
simple cycle of unitary operations, we can recover 
the familiar time-independent equations of contin- 
uum QM [61 [191 [301 El ig. Thus expUcit time de- 
pendence at the microscopic scale does not preclude 
compositions of unitary operations from being equiv- 
alent to more traditional time-independent models of 
physical dynamics. 

Even if there is nothing essential missing from ex- 
plicitly time-dependent QM models, it's not clear 
that they embody what we're looking for here: a 
classical-quantum correspondence that is so direct 
that we might hope to interpret the quantities and 
concepts of quantum dynamics in a simpler classical 
context. If we have to traverse a continuum limit 
in order to make contact with the familiar time- 
independent quantities and concepts of continuous 
dynamics, we lose this directness. 

1.4 Continuous dynamics 

Much of the simplicity of QM derives from the use 
of continuous operators. For example, if two QM op- 
erators A and B obey a commutation relationship of 
the form [A, B] = c for c a non-zero constant, then 
the operators A and B must necessarily be infinite- 
dimensionalj^ For this reason the fundamental oper- 
ators describing position and momentum in QM are 
continuous, just as they are in CM. 

When a (finite-state) physical system is described 
as a quantum system in continuous space, specifying 
the state of the system for a finite number of points 
in space completely determines the future evolution]^ 
This suggests that we may be able to construct con- 
tinuous CM models that are special cases of continu- 
ous QM dynamics if they also have this property. 

■^If j4 is a finite dimensional Hermitian operator with eigen- 
states \ai) and eigenvalues a^, and similarly for B, then c = 
{a,\AB - BA\a,) = (a; - ai)(a,\B\ai) = 0. 

finite-sized physical system with finite energy is effec- 
tively finite-dimensional, and the superposition of a finite num- 
ber of momentum eigenstates has only a finite number of in- 
dependent values in space |20j (see also Section |2.4|l. 



A famous example of a CM system of this 
sort is Fredkin's billard ball model of computation 
(BBM)[12]. This is an embedding of finite-state re- 
versible computation into continuous CM. At inte- 
ger times, all billiard balls are found at integer co- 
ordinates in the plane, all moving at the same speed 
in one of four directions. Collisions between billiard 
balls perform conditional logic — where a ball ends up 
after a collision depends on whether or not another 
ball was present to collide with itj^ 

1.5 Finitary classical mechanics 

A finite-sized BBM system is an example of what we 
might call a, finitary CM dynamics: one for which ex- 
act future states can be computed indefinitely using 
a finite number of logical operations. 

Finitary CM systems are special cases of CM that 
are equivalent to a finite state dynamics when ob- 
served only at integer times. In such systems the 
motion remains perfectly continuous between inte- 
ger times, but dynamical parameters (masses, spring 
constants, etc.) are chosen carefully and all dynam- 
ics is assumed to be perfect so that, as long as the 
initial state is chosen from a finite set of allowed 
states, one of the same allowed set is seen at each 
integer time (we only consider synchronous systems 
here, cf. [24) 

Since constraints are placed only on the values of 
dynamical parameters and the initial values of state 
variables, finitary CM systems map all of the con- 
tinuous Lagrangian machinery and conservations of 
CM onto finite-state dynamics At integer times the 

^This is an idealized model and is used here only to map fi- 
nite state computation into continuous CM and QM language. 
We are not concerned with issues of stability related to using 
only a set of measure zero of the possible initial states of a 
continuous system. 

^If, in the BBM, we want the allowed set of states to include 
all possible states with balls at integer positions moving at 
allowed velocities, then we need to modify the CM dynamics 
so that collisions (e.g., head on) that would take us off the 
integer grid do not occur — balls see zero potential in these 
cases and just pass straight through each other. 

^Least action principles involve continuous variations, and 
so there is not necessarily an equivalent principle that just in- 
volves the states at integer coordinates. The existence of such 
integer principles in single-speed systems (such as the ones 
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system is equivalent to a finite-state lattice dynam- 
ics such as a lattice gas automaton or partitioning 
automaton |32( I28j. From a continuum perspective, 
though, state variables arrive at integer coordinates 
in space at integer times — they don't remain on a 
fixed lattice at intermediate times. This is a subtle 
but important distinction. The same distinction in 
QM models allows infinite-dimensional position and 
momentum operators to be applied to particles that 
are always found at integer coordinates in continuous 
space at integer times. 




(a) 



(b) I • I 



1.6 Outline of the paper 

The remainder of this paper is divided into four sec- 
tions. 

Section Quantum concepts in a classical con- 
text discusses simple examples of recasting a classi- 
cal finite-state dynamics as a continuous finitary CM 
dynamics and thence as a continuous quantum dy- 
namics. The quantum description is closely related 
to classical sampling theory. Fermions and Bosons 
appear in the quantum description of classically iden- 
tical particles. 

5ec<«on[5[' Ideal energy uses the correspondence be- 
tween QM and finitary CM to define an ideal energy 
for finitary CM dynamics. This is the minimum to- 
tal energy that is compatible with the rate of local 
state change. A partitioning description of dynamics 
defines local state change. 

Section ^ An elastic string model and Section [5[' 
A colliding ball model define and analyze finite-state 
partitioning models that can be embedded in a con- 
tinuous relativistic dynamics. We discuss the rela- 
tionship between ideal energy and relativistic kinetic 
energy in these models. 



analysed in this paper) would not, however, contradict a re- 
cent no-go theorem |10 |. For models that have a macroscopic 
CM limit (such as the interacting examples we discuss), in- 
teger coordinate variations become infinitesimal variations in 
the limit. 



Figure 1: Embedding a finite-state classical dynam- 
ics, (a) Shift register with a single 1. (b) Classical 
particle in a box. (c) Quantum particle in a box. 

2 Quantum concepts 
in a classical context 

To illustrate some issues involved in describing CM 
dynamics using QM, and how QM concepts look in 
this classical context, we discuss in this section the 
simplest examples: non-interacting particle models. 
These models are embeddings of an N bit shift reg- 
ister. We start with the case where all bits of the 
shift register are O's except for a single 1 (Figure [T|i) . 
At each discrete time step the register is circularly 
shifted one position to the right: the 1 hops to the 
next position to the right, and if it is in the rightmost 
position it hops to the first position. 

2.1 Classical particle in a box 

The discrete-shift dynamics can be embedded into a 
continuous shift (Figure [T|d) . We replace the register 
with a one-dimensional space, and the 1 by a free par- 
ticle moving to the right at a constant speed s. The 
space is periodic, so that when the particle exits at 
one end, it reappears at the other end. If the length 
of the space is L, then we take our unit of distance 
to be L/N and our unit of time to be L/Ns. The 
particle is started at an integer coordinate, moving 
to the right at speed s. It will subsequently always 
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be found at one of the N integer coordinates at any 
integer time. Interpreting the presence or absence of 
a particle at a position as a 1 or a 0, the continu- 
ous system reproduces the discrete-shift dynamics at 
integer times. 

Conversely, the time evolution of the continuous 
system can be exactly predicted at integer times by 
evolving the discrete-shift system. This implies that 
the continuous system, with a finite set of allowed 
integer-time states, effectively has no more distinct 
states than the finite-state system. From the point 
of view of the discrete-shift dynamics intermediate 
states are fictitious states that have been added in 
order to construct a continuous embedding. 

In this example, as in the other finitary CM exam- 
ples we will consider, the finite subset of the contin- 
uous state set that can be visited by the dynamics 
at integer times is known when the system is defined 
and doesn't change with time. Periodic discreteness 
makes the details of the fictitious states between inte- 
ger times irrelevant for predicting the long-term be- 
havior of the dynamics. 

2.2 Particle in a quantum box 

Given finite energy a QM particle in a box can be 
completely described using only a finite set of dis- 
tinct states. One could regard the CM embedding 
of the shift register example above as a semiclassical 
analog of a QM particle in a box, illustrating how 
a system with a dynamics that obeys a continuous 
translational symmetry can be made to act like a 
finite state system by restricting the set of allowed 
initial states in a manner that is conserved at integer 
times by the dynamics. 

We can make this more than an analogy by recast- 
ing the continuous CM "particle in a box" dynamics 
in QM language. 

2.2.1 Unidirectional Hamiltonian 

Consider a free massless relativistic particle in a one 
dimensional periodic box. If we ignore for a moment 
that the particle can travel both left and right we 
can take the Hamiltonian to he H = cp, where p = 
^ihd/dx is the usual momentum operator and c is 



the speed of light. The Schrodinger equation becomes 

(1) 



ih— = Hip = -cih— 
ot ox 



This has the solution ^{^{x, t) = f(x — ct), where f{x) 
is the initial state of at i = and can be specified 
arbitrarilyj^ This dynamics just shifts whatever state 
we start with uniformly to the right at speed c. 

Now we assume that this system has finite energy. 
In QM, this assumption imposes constraints on the 
allowed initial states. The conventional way to ana- 
lyze this for a QM particle in a box is to describe the 
particle's position in terms of momentum eigenstates: 
only waves that fit an integer number of wavelengths 
across the width of the box are allowed by the bound- 
ary conditions. Since momentum eigenstates with 
higher spatial frequencies correspond to higher ener- 
gies, a finite maximum energy eigenstate only allows 
us to use a finite number of frequency components. 

2.2.2 Energy based counting of states 

The QM average energy corresponds to the classi- 
cal energy, and so the natural constraint for a closed 
classical dynamical system is that the average energy 
is bounded (rather than the maximum energy eigen- 
value) . We will take the average energy constraint as 
the starting point for the analysis in this section. 

Average energy is related to rate of state 
change [27]: for any closed QM system, the average 
energy E (taking the ground state energy as zero) 
gives an achievable bound on the rat^ i'± at which 
the system can pass through a long sequenc^^of dis- 
tinct (i.e., mutually orthogonal) states: 



v± < 2E/h 



(2) 



*The close connection between classical and quantum for 
this Hamiltonian is discussed in the context of hidden variables 
theories in |37| . The embedding of classical EM into a QM 
description in [7] is also closely related. 

^It is easy to see that orthogonality (or approximate orthog- 
onality) must always occur at regular intervals, since the inner 
product of any two states \ipt) = Ut\^o) and I'/'t+s) = Utl^s) 
separated in time by an interval s depends only on s. 

^'^The bound actually depends slightly on the length n of 
the total-system dynamical cycle that an orthogonal change is 
part of — we omit a factor of n/ (n — 1) on the right side of the 
inequality. 
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Now suppose the particle is initially very well local- 
ized. Then the initial wavefunction f{x) is sharply 
peaked at one point and zero everywhere else. Un- 
der the "shift" dynamics, after a short time the sharp 
peak will have moved enough so that the old and new 
position states will be distinct (Figure . The more 
sharply the particle is localized, the sooner the states 
will be distinct and so the higher the rate of orthogo- 
nal change and hence the higher the average energy. 
If instead the "peak" is spread evenly over the whole 
box, the shift will never produce a different state and 
the average energy can be very low. 

Given a finite average energy E there is thus a 
limit to how sharply the particle can be localized and 
this determines how many distinct position states are 
available to the system. If the width of the box is L 
and there are N distinct position states, then the sys- 
tem transitions to a new orthogonal state in the time 
it takes the pattern to shift a distance of L/N. Thus 
the rate at which the system passes through orthogo- 
nal states is = cN/ L and so, from Equation |2] N 
can't be more than 2LE/ch. The maximum number 
of distinct states is bounded by the energy E. Since 
the average (classical) momentum is E/c, this can 
also be written as < 2pL/h. 



2.2.3 Bidirectional Hamiltonian 

The unidirectional Hamiltonian we've used here 
seems unphysical, since a real particle should be able 
to move in either direction. The direction should be 
part of the state, not part of the Hamiltonian. This 
can easily be fixed. Since the Hamiltonian H — cp 
moves the particle uniformly to the right at speed 
c without changing the shape of the wavefunction, 
clearly the Hamiltonian H = —cp would carry the 
particle uniformly to the left at speed c. If we make 
a two component state vector and use the Hamilto- 



H = 



cp 






-cp 



(3) 



represent it in the second component If we define 

(4) 



1 

-1 



then the Hamiltonian can be written as H — cap. 
Thus (T can be interpreted as the operator that reads 
the sign of the particle direction from the state infor- 
mation. This is a one-dimensional version of the Weyl 
equation, sometimes also referred to as the massless 
Dirac equationj^ 

2.3 Discrete position basis 

In Figure [TJ;, we depicted the wavepacket represent- 
ing the particle as a square wave. In this section we 
discuss the actual shape of the wavepacket needed 



then the particle will travel to the right if we repre- 
sent it in the first component, and to the left if we 



to achieve the energy bound of Section 2.2.2 on the 
number of distinct states. 

If we were just specifying a discrete distribution at 
integer times and positions, the probabilities needed 
to reproduce the shift dynamics are obvious: the par- 
ticle is at some integer position with probability 1 
and at other integer positions with probability 0. We 
want to extend this distribution to the continuum. 

For continuous positions, the width of the 
wavepacket represents uncertainty in where exactly 
the particle would be found if we looked. Since the 
shape of the wavepacket doesn't change as it shifts, 
we are free to interpret all of the uncertainty as un- 
certainty in the initial position of the particle. 

There is, however, a better way to interpret this. 
The shape of the wavepacket is chosen to allow us 
to represent exactly N distinct states, each with the 
wavepacket centered at an integer coordinate. The 
uncertainty in particle position blurs distinctions, 
combining continuum states that would otherwise be 
distinct. Thus another way to interpret the posi- 
tion uncertainty would be to say that we invented 

-'^-'^The energy associated with this Hamiltonian is always pos- 
itive as long as we use only positive momentum eigenvalues to 
describe motion to the right at speed c, and only negative mo- 
mentum eigenvalues for motion to the left. 

^^Lattice models that reproduce the Weyl or Dirac equation 
in the continuum limit are discussed in [6l ll9l[30ll9ll43) . Note 
that here we do the opposite: we use the Weyl equation to 
exactly simulate discrete lattice dynamics at integer times. 
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(a) (b) 

Figure 2: Cardinal Sine Probability Distribution, (a) Normalized sinc^ distribution function (solid) versus 
normalized Gaussian of the same height (dotted), (b) Real (solid) and imaginary (dotted) components of 
sinc-based amplitude distribution. 



fictitious intermediate positions in order to map a 
finite-state dynamics onto the continuum, and then 
we treated the exact value of the position coordinate 
as if it was uncertain in order to eliminate the ficti- 
tious position information by representing it proba- 
bilistically as missing information. 

Now, if we represented the particle by a square 
pulse, as in Figure [T]:, we would be assigning equal 
probabilities to all fictitious intermediate positions in 
one region and zero elsewhere, with a sharp transi- 
tion in between. The only constraint on the con- 
tinuous wave function, however, is that it should ex- 
actly represent the finite-state dynamics when sam- 
pled at integer times and positions. We should there- 
fore choose the wavefunction that includes as little 
information as possible about fictitious intermediate 
positions (i.e., in keeping with Jaynes Maximum En- 
tropy Principle |17j . we should choose the wavefunc- 
tion that is maximally noncommital with regard to 
missing information). This implies a bound on the 
Fourier spectrum of the wavefunction: spatial fre- 
quencies that are needed only to specify details about 
fictitious intermediate positions should be omitted. 

In Figure we illustrate a probability distribution 
function based on the cardinal sine function, 

sin a; , , 

smc X = (5) 

X 

Cardinal sine gets its name as the basis of E. T. Whit- 
taker's cardinal interpolation function, and it plays a 



central role in bandlimited sampling theory |35l 131] . 
It is the Fourier transform of a finite range of fre- 
quency components, all taken with equal weight. The 
probability distribution function shown with the solid 
line is sinc'^Trit. This is a normalized distribution that 
is close to a Gaussiarp^ (dotted line) near the mean, 
has the value 1 at the point u = and has the value 
at every other integer position. 

The corresponding amplitude distribution for 
a particle centered at integer position m is 
sine 7r(M — m) times a phase factor. Of course this 
function also has magnitude one at the center of the 
distribution, and magnitude zero at all other inte- 
ger positions. We can use a weighted sum of such 
functions to construct any assignment of amplitudes 
to integer positions. Figure [2]d shows the real and 
imaginary parts of the function 

(/)(7rM) = e'''" sine nu (6) 

The set of functions (f)(Tr(u — m)) for m an integer 
are all mutually orthogonal: these functions form a 
complete basis. We will refer to these states as the 
discrete position basis. 

The discrete position basis states (and any states 
constructed out of them) have a bounded spatial fre- 
quency spectrum. These basis states answer the ques- 

^''The use of low-pass filtering to extend a distribution with 
sharply localized features into a smooth continuous one is in- 
timately related to diffusion |21) . 
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tion we started this section with: for the shift dy- 
namics, they have the waveshape that achieves the 
maximum possible rate of orthogonal change for a 
given average energy (Equation [2|. The average en- 
ergy of a particle in a discrete position basis state 
is /i/2r, where t is the time it takes the particle to 
travel from one integer position to the next. We de- 
rive these fastest-changing states directly from the 
Hamiltonian H = cp in the next section. 

2.4 Bandlimited dynamics 

A general construction for fastest changing states is 
given in p7l. These states are essentially Fourier 
transforms of a finite set of energy eigenstates, which 
are the slowest changing states. 

In the macroscopic limit, the discrete position basis 
states of Equation [6] are fastest changing states (for 
given average energy) of a ID particle evolving under 
the shift dynamics. The case we are interested in, 
however, is finite and periodic. It is instructive to 
construct the exact states that achieve the energy 
bound of Equation [2] in this case. If we consider a 
situation where the particle travels only to the right 
(Equation [ij , then the energy eigenstates satisfy the 
equation 

-'"^ (7) 



cih 



dx 



We assume that our energy scale is chosen so that 
the smallest value of E is zero. This equation is sat- 
isfied by an exponential of x. If we use positive x 
coordinates that range from to i in the box, then 
the generic solution that is periodic at a; = i is given 
by il^n oc e^'^™^/^, where n is a non-negative integer. 
This function satisfies Equation [7] provided 



E = nh- 



L 



(8) 



In our CM embedding of an iV-bit shift register 
dynamics, the only meaningful positions at integer 
times are x^ — uL/N for u an integer: the wave- 
function at the (fictitious) non-integer positions is a 
continuous extension, completely determined by the 
values at the integer positions. Thus we require that 
the probabilities of finding the particle at integer po- 



sitions at integer times add up to one, and so 

1 



■>Pn{Xu) 



N 



(9) 



for which '^^il'ni^u)4'n{xu) — 1- Using Equations [S] 
and [9] we can construct a normalized fastest changing 
state centered at position zero: in this case it is sim- 
ply the sum of the N different ipnS taken with equal 
amplitudes [22]. Therefore the fastest changing state 
centered at position is given by 



V ^ ,1=0 
N-1 

_ _j_ \ " 2TTin{u~m)/N 

- N ^ 

n=0 



(10) 



(11) 



sin7r(M-m) 



(12) 



where the final expression comes from summing 
the geometric series and some rearrangement. The 
fastest changing states are Fourier transforms of the 
energy eigenstates. For large TV, ^mixu) turns into 
cj){Tr{u—m)) of Equation|6] <&m(a:«) also has the prop- 
erty that it is one for u = m (mod N) and zero for u 
any other integer. For integer values of u, the ^m{xu) 
functions provide a complete discrete-position basis 
for the finite system. 

2.5 Both continuous and discrete 

The discrete position basis state ^m{xu) corresponds 
to a probability distribution that is normalized both 
from a continuous and a discrete point of view. Fur- 
thermore, the discrete normalization doesn't apply 
just to a sampling of the state at integer locations: 
it applies to any discrete sampling of the state at a 
unit-spaced set of positions. 

The continuous normalization is apparent from 
Equation |11[ since 

dx %{x) <i>o{x) 



^ N-1 N-1 



iV2 



2. 



n=0 n'=0 



N 



(13) 
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To see the generality of the discrete normalization, 
we similarly sum over positions Xu+s where u is an 
integer and < (5 < 1. This shifted normalization 
also implies that if we start the system in a discrete- 
position basis state, then as the wavefunction shifts 
continuously with time the values sampled at fixed 
integer locations always constitute a normalized dis- 
crete representation of the state. 

Both normalization properties are of course in- 
herited by any normalized wavefunction \['(a;,t) con- 
structed out of the states (i.e., any \E' whose spec- 
trum of spatial frequencies is bandlimited like the 
<&'s). 'i'{x,t) can be exactly reconstructed from the 
amplitudes ^(xfe,t) sampled at integer locations: 



(14) 



fc=0 



2.6 Intermediate positions 

We can construct a <i>c state centered at any contin- 
uum position c as a superposition of the N discrete- 
position basis states centered at integer locations 
m. This is obvious since the states form a com- 
plete basis for describing the continuously shifting 
time evolution of an initial state. A shift of all 
states by the same amount yields a new basis. 
If we start our system with a particle at an inte- 
ger location in a basis state and then change bases 
in step with the particle motion, we can follow the 
continuous motion of the particle, always seeing it in 
a single localized basis state. From this point of view 
the quantum motion is just as spatially local and con- 
tinuous as the classical motion: in both cases we see 
the constraint on the classical system of a discrete 
set of possible positions at integer times reflected as 
a corresponding discrete set of possible intermediate 
positions at each intermediate time. 

2.7 Classical Fermions and Bosons 

The shift register example is easily extended to mul- 
tiple particles. This simply corresponds to a value in 
the register (Figure [T^) with more than a single '1' 
bit. For simplicity, we'll assume for the moment that 
all particles are traveling to the right. 



In elementary QM, the overall state of a collection 
of independent particles is normally represented us- 
ing a product of single-particle states: we track the 
state of each particle separately as if it were the only 
particle in the system. 

For our classical shift register we've already de- 
scribed a single particle system using a localized 
wavepacket basis (Equation 12 1. Here we'll denote 



the single particle basis state where particle j has its 
wavepacket centered at location x as \x,j). An over- 
all state can be described as a product of such states, 
one per particle. 

Since all particles are moving to the right, the di- 
rection CT is -1-1 in all cases. The Hamiltonian for 
a multi-particle system is then H — cpj , where 
Pj acts only on particle j. For example, for a two 
particle system (i.e., shift register with two I's), the 
Schrodinger equation would be 



ih 



d\x,0)\y,l) 



dt 



{cpo + cpi)\x,0)\y, 1) 
(cpo|x,0))(cpi|y,l)) (15) 



This equation is satisfied if both ihd\x,0)/dt = 
cpo\x,0) and ihd\y,l) /dt — cpi\y,l) hold, and so the 
two particles each follow a shift dynamics. 

By using a product of single-particle wavefunctions 
to construct our multiparticle wavefunction we've 
added two kinds of extra states that didn't exist in 
the original shift register. First of all, our states can 
represent more than a single 1 at the same location, 
which isn't possible in the original shift register. Sec- 
ondly, our states keep track of which particle is which, 
but these particles are identical I's in a shift register 
and interchanging two I's in a binary number doesn't 
give us a different state! 

We can fix both of these problems by combining 
each set of equivalent states into a single state. We do 
this by adding together all product states that differ 
only by a relabeling of identical particles and we call 
that a single state. In performing this sum on a set of 
M equivalent states we weight half of the states by a 
factor of +1/\/M and half by -1/\/M. This is done 
in such a manner that the entire combination changes 
sign if we interchange any two particle labels (i.e., it 
is antisymmetrized) . In our two particle example, the 
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antisymmctrized wavefunction is 

|x,0)|y,l)-|y,0)|.T,l) ^^g^ 

Using only antisymmctrized wavefunctions we can no 
longer represent two particles (I's) at the same loca- 
tion: as in this example, if two mean positions are 
equal (i.e., x = y) the wavefunction is zero. 

Replacing each set of equivalent states with a single 
antisymmctrized sum avoids over-representing states. 
For a classical embedding into QM, a superposition 
of classically distinct states corresponds to a classi- 
cal statistical ensemble: we can interpret all prob- 
abilities as arising from ignorance about which sin- 
gle state the system started in. Thus by assigning 
equal probabilities to a set of equivalent states, the 
antisymmctrized superposition represents the equiv- 
alent states as if we were completely ignorant about 
which particle was which. Of course the particle label 
information doesn't exist in the original finite-state 
dynamics and so, once again, the probabilities we've 
assigned represent ignorance of fictitious classical de- 
tails added in constructing our QM description. 

Thus for a classical dynamics of bits, if a fixed num- 
ber of identical I's are described quantum mechani- 
cally as if we could distinguish which is which, and 
then the fictitious distinctions are eliminated using 
superposition, they are Fermzonsj^ If, instead, there 
were no constraint in our original classical dynamics 
on how many I's can be represented at each location, 
we would use symmetrization rather than antisym- 
metrization to eliminate the fictitious distinctions. In 
this case the I's would be Bosons. 

2.8 Non-classical amplitudes 

There are two kinds of probabilities that can arise in a 
quantum description of a classical finite-state dynam- 
ics. There may be unknown information about the 
state of the actual dynamics: this is represented by 
ordinary probabilities and contributes to the entropy 
of the system. There may also, however, be unknown 
information about fictitious classical details that were 

'^^Deterministic Fermions have been described previously in 
a formal QM field-theoretic context (e.g., |36) V 



added in constructing the quantum description of the 
dynamics (such as distinct labels for identical I's): 
this is represented by probabilities derived from am- 
plitudes and doesn't contribute to the entropy of the 
system. 

There is more to amplitudes than this, though. If 
the dynamics was always described in a single clas- 
sical basis, then the superposition would evolve like 
a classical statistical ensemble and the coefficient of 
a basis state in the superposition could be any func- 
tion of the basis state's ensemble probability — there 
is no recombining of coefficients in the ensemble dy- 
namics. The fact that amplitudes are square roots 
of probabilities makes it possible to also use them 
to describe the system in multiple bases: the sum of 
the squares of component magnitudes is the same for 
a vector described in any basis. Amplitudes, when 
used for basis change, still don't represent unknown 
information about the real system, and so again don't 
contribute to its entropy. 

2.9 Finite-state unitary dynamics 

For a classical finite-state system with a conserved 
number of identical I's, we can treat the I's as dis- 
tinguishable particles in a QM description and follow 
their trajectories. We can then use appropriate sym- 
metrization of superpositions of products of single 
particle states to merge equivalent states that aren't 
actually distinguishable. 

We can extend this approach to describe classi- 
cal systems in which the numbers of I's and O's are 
not constant. The creation (and annihilation) oper- 
ators of quantum field theory allow us to act on an 
appropriately symmetrized superposition of products 
of single-particle states, to add (or remove) a single- 
particle basis state in each product term while main- 
taining the symmetrization of the overall state. We 
can use an occupation-number basis to avoid deal- 
ing directly with single-particle states, but a vestige 
of antisymmetrization persists as anticommutation of 
Fermionic field operators |44|. 

We can avoid all of this complexity entirely, 
though, by simply not describing identical I's as if 
they were distinguishable! We won't have to nullify 
fictitious particle labels using symmetrization if we 
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don't put them in in the first place. Similarly, unnec- 
essary complexity is added when we employ a con- 
tinuum of positions and times to describe a classical 
system with only a finite number of possible states. 

Both sources of complexity are eliminated if the 
QM dynamics is constructed to be isomorphic to the 
classical finite-state dynamics. If we start with a clas- 
sical dynamics expressed as a composition of invert- 
ible logic operations we can, in the QM description, 
simply replace each classical logic operation with an 
equivalent finite-time finite-state unitary transforma- 
tion. There are then no extra states added to the de- 
scription that need to be nullified. Moreover, if the 
finite-state dynamics can be interpreted as an integer- 
time sampling of a finitary CM dynamics, then it is 
easy to identify the familiar continuous-space and - 
time conserved quantities of CM, including classical 
energy. The corresponding identification of QM en- 
ergy is discussed in the next section. 

3 Ideal energy 

In CM there is no minimum energy required to allow 
a given sequence of state changes to take place in a 
given period of time. In QM, there is. If a classi- 
cal finite-state dynamics is recast isomorphically as 
a unitary QM dynamics we can place a bound on 
the least possible average energy for physically real- 
izing the classical dynamics — the ideal energy of the 
dynamics (cf. [21]). 

For a dynamics described as a composition of 
spatially-localized invertible operations. Equation [2] 
allows us to define a minimum energy bound 
£'statc-change that depends on counting state changes: 

£^statc-chango = * (state-chaugcs in r) (17) 

ZT 

Here t is the length of the time interval during which 
we count state changes. Each invertible operation 
that actually changes the state counts as one state 
change. 

Now, in CM, energy of motion appears as local- 
ized state change. Thus the state-change energy we 
define here can be thought of as a kind of kinetic en- 
ergy: it bounds classical energy of motion, at a given 



time. The total classical energy is then bounded by 
the maximum value that the state-change energy ac- 
quires in the course of the time evolution. 

In this section we discuss state-change energy and 
use it to define the ideal minimum energy of a finitary 
CM dynamics. In Sections |4] and |5]wc apply this idea 
to two examples of finitary CM dynamics. 

3.1 Classical partitioning dynamics 

If a spatially extended finite-state system evolves in 
time synchronously, the rate of global state change 
dramatically underestimates the minimum energy 
needed for an actual physical realization of the dy- 
namics. In particular, the energy of a collection of 
independent subsystems is additive: localized inde- 
pendent changes must be counted separately. 

We can arrive at a more realistic minimum-energy 
estimate by making use of the locality and time de- 
pendence inherent in a physical description of clas- 
sical finite state dynamics. Consider, for example, 
a finitary CM dynamics such as Fredkin's BBM. In 
this model, all collisions happen synchronously, at a 
discrete set of possible locations at a discrete set of 
possible times. If we focus on a time interval start- 
ing shortly before a possible collision time and ending 
shortly after, then during this interval the system can 
be partitioned into a set of disjoint regions, each of 
which evolves independently: within each region one 
or more particles head towards the locus of a possible 
collision, collide or don't collide, and then head away 
from the locus. Which regions evolve independently 
must of course change with time or the particles could 
never cross between regions. Thus the partitioning 
used in a description of the system as a collection 
of independent subsystems must be time dependent, 
even though the continuous dynamics itself follows a 
time independent dynamical law. 

Figure [3] illustrates the partitioning idea applied 
to the "classical particles in a periodic box" model 
discussed in Section |2] A time history from this fini- 
tary CM dynamics is shown, with the positions of two 
continuously moving particles (I's) sampled at inte- 
ger times. If the spatial intervals marked with solid 
lines in the figure are taken to be blocks of a partition, 
then during the period between one integer time and 
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of one partition, and then on each block of the other 
partition. 



Figure 3: Integer-time history of two particles moving 
continuously in a periodic box. 



the next any particles that start off in a block remain 
in that block and only interact with other particles 
in the same block. The classical finite-state dynam- 
ics at integer times can be summarized as: swap the 
values (I's or O's) at the two integer locations in each 
block of the partition. 

Notice that, looking just at the integer-time his- 
tory, it is ambiguous whether continuously moving 
particles whose paths intersect pass through each 
other or bounce back. This kind of ambiguity is very 
common in extending a finite-state dynamics to con- 
struct an equivalent finitary CM dynamics. If, how- 
ever, we only want to quantify the minimum possi- 
ble amount of motion consistent with a given classi- 
cal finite-state dynamics, then there is no ambiguity. 
For example, for the dynamics illustrated in Figure[3] 
the least amount of motion occurs if the two particles 
don't move between steps 2 and 3. 

Of course in a real CM dynamics, the amount of 
motion that occurs inside the block where the two 
particles collide and bounce back could be small if 
they slow down quickly, but not exactly zero. We 
might think of zero motion either as a limit or, alter- 
natively, as the exact amount of motion seen at the 
midpoint of the block update. 

3.2 Ideal kinetic energy 

The integer-time swap dynamics of Figure [3] can be 
described as a time-dependent sequence of unitary 
operations, with operations acting first on each block 



Equation 17 assigns h/2T of state-change energy 
to each block update involving a single moving par- 
ticle. This agrees with our earlier QM analysis of the 
minimum achievable energy for a freely moving CM 
particle. Since there is no kinetic energy associated 
with blocks that don't change, either in the CM or 
the QM description, the state-change energy defines 
an ideal kinetic energy in this case: it is at all times 
the minimum kinetic energy that could possibly be 
achieved in a physical realization of the dynamics. 
Its maximum value is the ideal total energy. 

In general, the definition of kinetic energy in a fini- 
tary CM dynamics is somewhat ambiguous. For ex- 
ample, the same dynamics might be interpretable as 
either a relativistic system or a non-relativistic sys- 
tem 121]. If the state-change energy is taken as the 
ideal kinetic energy at all times, this ambiguity is 
resolved. Even if this is not done, however, the state- 
change energy sets the ideal energy scale, since it 
bounds the minimum total energy that can be in- 
volved in each possible block transition: for any CM 
interpretation of the dynamics one of these bounds 
will be the most constraining. 

Below, we discuss ideal energy in two example 
models, beginning with a model that is closely re- 
lated to the swap dynamics. 



4 An elastic string model 

In this section we discuss a classical finite state model 
of wave motion in an elastic string. The finite-state 
string model corresponds exactly to the integer-time 
behavior of a continuous CM model of the type dis- 
cussed in Section |1.5[ The finite-state string model 
also exactly reproduces the behavior of the time- 
independent one-dimensional wave equation sampled 
at integer-times and locations. There is a finite-state 
constraint on the allowed initial shape and motion 
of each unit segment of the wave and this constraint 
reappears at each integer time. In the continuum 
limit the discrete constraint on the wave disappears; 
the exactness of the wave dynamics itself (at discrete 
times) is not dependent on this limit. 
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Figure 4: Discrete wave dynamics. Elastic string is held fixed where it crosses black bars. 



The model presented here has been discussed be- 
fore [ini HSl US] , but the analysis of translational mo- 
tion, the relativistic interpretation and the analysis 
of ideal-energy given here are all new. 

4.1 Discrete wave model 

Consider an ideal continuum string for which trans- 
verse displacements exactly obey the wave equation. 
In Figure |4|l we've illustrated an initial configuration 
with this string stretched between equally spaced ver- 
tical bars. The set of initial configurations we're al- 
lowing are periodic, so the two endpoints must be 
at the same height. Any configuration is allowed as 
long as each segment running between vertical bars 
is straight and lies at an angle of 45°. 

Initially the string is attached at a fixed position 
wherever it crosses a vertical bar. We start the dy- 
namics by releasing the attachment constraint at all 
of the gray bars. The attachment to the black bars re- 
mains fixed. In Figure|4]3 the segments that are about 
to move are shown with dotted lines: the straight seg- 
ments have no tendency to move. Under continuum 
wave dynamics, the dotted segments all invert after 
some time interval r. This will be our unit of time 
for the discrete dynamics. The new configuration at 
the end of this interval is shown in Figure |4|;. At 
this instant in time all points of the string are again 
at rest and we are again in an allowed initial config- 
uration. At this instant we interchange the roles of 
the black and gray bars and allow the segments be- 
tween adjacent gray bars to move for a time inteval 



r. The dynamics proceeds like this, interchanging the 
roles of the black and gray bars after each interval of 
length T. Since attachments are always changed at 
instants when all energy is potential and the string 
is not moving, the explicit time dependence of the 
system doesn't affect energy conservation. 

We express this dynamics as a purely digital rule 
in Figure [Sj In Figure [5^ we show the wave just af- 
ter the evolution from the previous figure, with the 
black bars shown marking the attachments for the 
next step. To simplify the figure we have suppressed 
the gray bars — they are always situated midway be- 
tween the black bars and so don't need to be shown. 
We have also added a grid of 45° dotted lines that 
shows all of the segments that the string could pos- 
sibly follow. In Figure |5]3 we add in horizontal black 
bars, in order to partition the 2D grid into a set of 
2x2 blocks that can be updated independently. Note 
that in all cases the segments that are allowed to 
change during this update step, as well as the cells 
that they will occupy after the update, are enclosed in 
a single block. Segments that aren't going to change 
stretch across multiple blocks. The ID box below 
Figure [5]d contains just the slope information from 
the string. This array of gradients is clearly suffi- 
cient to recreate the wave pattern if the height at 
one position is known. 

Figure [5]: shows the dynamical rule for a block. 
Since the dotted lines indicate the direction in which 
segments must run if they appear in any cell, the 
state information for each segment is only whether 
it is there or not: this is indicated with a or a 
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Figure 5: Discrete wave dynamics, (a) Wave configuration with possible wave paths indicated (dotted lines), 
(b) One of two partitions used for a discrete update rule (horizontal and vertical lines). Array of wave 
gradients shown at bottom, (c) Top, dynamical rule for wave. Presence of wave-path segments is indicated 
by I's. Bottom, equivalent dynamical rule for gradients. 



1. The only segments that change are peaks /\ or 
valleys \/, and these are represented by two I's at 
the top of a block or at the bottom of a block re- 
spectively. The rule is that peaks and valleys turn 
into each other, and nothing else changes. This rule 
applies to the blocking shown, and then to the com- 
plementary blocking when the attachments change, 
alternating back and forth. This is a partitioning dy- 
namics. 

4.2 Exact wave behavior 

At the bottom of Figure |5]: we've presented a dynam- 
ics for the gradient of the wave. Since the full 2D dy- 
namics just turns peaks into valleys and vice versa, 
leaving straight segments unchanged, we can perform 
that dynamics equally well on the array of gradients! 
As the 2D dynamics interchanges which blocking to 
use, the dynamics on the gradients also alternates 
which pairs of gradients to update together. In all 
cases, the dynamics on the gradients duplicates what 
happens on the string: if the two dynamics are both 
performed in parallel, the gradient listed below a col- 
umn will always match the slope of the string in that 
column. 

The dynamics on the gradients has a very inter- 
esting property. Turning a peak into a valley and 



vice versa is exactly the same as swapping the left 
and right elements of a block. Leaving a // or \\ 
unchanged is also exactly the same as swapping the 
left and right elements of a block. In all cases, the 
dynamics on the gradients is equivalent to a swap. 

This means that the left element of a block will 
get swapped into the right position, and at the next 
update it will be the left element of a new block and 
will get swapped into the right position, and so on. 
Thus all of the gradients that start off in the left side 
of a block will travel uniformly to the right, and all 
of the gradients that start in the right side of a block 
will travel uniformly to the left. 

This shows that the system obeys a discrete version 
of the wave equation. Half of the gradients constitute 
a right-going wave, and half constitute a left-going 
wave. At any step of the dynamics, the 2D wave 
in the original dynamics is just the sum of the two 
waves: it is reproduced by laying gradients end to 
end! 

As the number of cells in our lattice gets arbi- 
trarily large, the straight segments in our waves be- 
come smaller and smaller compared to the total width 
of the lattice, and the discrete wave equation turns 
smoothly into exactly the continuum wave equation. 
As we will see, though, even the discrete model ex- 
actly obeys the continuum wave equation. 
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Figure 6: Rescaling a discrete wave, (a) Original wave and array listing its gradients, (b) Each pair of 
gradients is repeated four times. The corresponding wave is drawn by placing gradients end-to-end. (c) 
Equivalent superposition of continuous waves. 



4.3 Rescaling limit 

The discrete wave system also has another contin- 
uum limit that allows us to embed our discrete wave 
dynamics with time-dependent blocking into a con- 
tinuous dynamics with a time-independent dynamical 
law. 

Take an array of gradients that define a wave (Fig- 
ure [6^) and make a longer list by replicating every 
pair of segments N times, producing a list that is N 
times as long (Figure ^jp). Since the dynamics for 
the gradients is simply that all of the even numbered 
gradients shift one way and all of the odd numbered 
gradients shift the other, the state after N steps cor- 
responds exactly to the state we would have gotten by 
running one step on the original configuration, pair- 
ing the gradients into the alternate set of blocks and 
then replicating each block N times. 

Although the iV-step dynamics is exactly equiva- 
lent to the original, the waves we reconstruct from 
the array of repeated blocks look much different than 
simply a scaled up version of the original wave: a 
peak /\ doesn't turn into a peak times as large, 
it turns into N little peaks /\/\/\ . . ./\. This is 
a wiggly flat line, and as N grows the wiggles get 
proportionately smaller and smaller compared to the 
straight parts // and \\, which simply get N times 
longer. In the limit, the peaks turn into flat places 
with a net velocity up or down, still obeying the con- 



tinuum wave equation. 

4.4 Continuous at all scales 

The rescaled continuum configurations correspond 
isomorphically with the discrete time evolution, since 
a fiat portion of the wave moving up or down can 
simply be interpreted as an array of peaks or val- 
leys. In fact, if we simply draw each block of the 
active partition that we've been showing as /\ or \/ 
as a horizontal segment moving up or down instead, 
then the discrete dynamics is directly mapped onto 
the continuous wave equation without invoking any 
limiting process! 

Since this evolution obeys the continuous wave 
equation we can draw it as the superposition of con- 
tinuous right going and left going waves. This is shown 
in Figure The top solid wave is a redrawing of 
Figure [6^ with a horizontal segment used instead of 
the /\ (or equivalently, the continuum limit of Fig- 
ure ^jp) . It is also the sum of the dotted rightgoing 
and left going waves shown below it, each segment of 
which has slope ±1/2. The collection of gradients 
that make up the rightgoing and leftgoing waves are 
shown at the bottom. These are just the even posi- 
tion (top) and odd position (bottom) gradients from 
Figure [6|i, each stretched to fill the width of a block. 
After the waves each move half the width of a gradi- 
ent segment, all segments will again be aligned and 



15 



will add up to a result that corresponds to the next 
step of the discrete wave dynamics, with the appro- 
priate blocking. 

4.5 Relativistic motion 

Assume the string carrying a discrete wave wraps 
around a space of width N. We've discussed the 
horizontal motion of waves along such a string, but 
the string itself can move vertically. For example, a 
pattern such as /\/\/\. . ./\ all the way around the 
space reproduces itself after two steps, but shifted 
vertically by two lattice units. This is clearly the 
maximum rate of travel for a string: one position 
vertically per time step. 

We can express the net velocity of the string in 
terms of the populations of rightgoing and leftgoing 
gradient segments. With N positions there are N/2 
rightgoing segments and N/2 leftgoing segments, and 
the members of each group never change. Thus if 
there are R rightgoing \'s, there must be {N/2) — R 
rightgoing /'s, and similarly for the other direction. 
Therefore if we know the numbers R and L of right- 
going and leftgoing \'s, we have complete population 
information. 

For a wave that wraps around the space and joins 
at the ends, there must be equal numbers of /'s and 
\'s. Thus half of the gradient segments must be \'s. 



R + L = 



N 



(18) 



The average upward velocity of the string depends on 
the difference between R and i. 



R- L 
N/2 



(19) 



If we think of the dynamics as swapping adjacent 
pairs of gradients during each step (even though \\ 
and // pairs don't actually change), each gradient 
moves one position up or one position down during 
each update. A \ moving right always moves the 
gradient it passes one up as it goes by, whereas a \ 
moving left moves the gradient it passes one down. 
Thus we get Equation [l9] for the average velocity of 
gradients that pass the N/2 \'s. The average veloc- 
ity of gradients that pass the /'s is exactly the same 



(their populations are complementary but the effects 
of passing them are also complementary) . 

Since the motion of each gradient contributes a 
velocity of ±1 to the the net motion of the string, 
it is natural to interpret the net motion relativisti- 
cally. We can focus on just the \ particles, since 
the / particles have the same behavior. The to- 
tal energy is proportional to the number of identi- 
cal particles all moving at the same speed, so we let 
E ^ R+ L 
(taking c 



N/2 (Equation p|. Since P = Ev 
1), this gives us a vertical momentum 



oi P = R — L (Equation 19 1. Then the mass of the 
string must be AP = E^ - P^ = 4RL and E = 7M 
where 7 = (1 — v^)^^/^. 

We can exactly analyze the vertical motion of the 
string statistically, since the probability p — R/{R + 
L) is the exact frequency with which leftgoing gradi- 
ents will encounter a rightgoing \ in the course of a 
cycle of length N/2. 

A rightgoing gradient is the left element of a block 
of two adjacent string segments that are about to be 
updated by the dynamical rule. If the block contains 
\/ then the dynamics will move both elements up- 
ward one position. Thus p is the probability that, for 
any block holding part of the string, the left element 
of the string is ready to take an upward step, and 
I — p the probability it is ready to take a downward 
step. 

Similarly, if a block contains \ as its right element 
(i.e., a leftgoing \), then the block is ready to take a 
downward step if the left element is /. The chance 
that the block contains a leftgoing \ is L / {R + L) 
which is just 1 — p. Thus in either position of a block, 
p is the probability a gradient is ready to take an 
upward step and 1 — p the probability it is ready to 
take a downward step, and so the average frequency 
(per step) of both elements of a block moving upward 
is and both downward is {1 — p)^. The average 
upward velocity of the string is thus — (1 — p)^ 



which just gives us back Equation 19^ 



^^Since p is also the probability of moving a gradient upward 
as it passes and 1 — p the probability of moving it downward, 
the average vertical velocity v is also given hy v = p — (1 — p). 
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4.6 Relativistic kinetic energy 



4.7 Time dilation 



Our original continuous description of the elastic 
string dynamics involved minimal- motion between in- 
teger times: all cases remain motionless except for 
\/ and /\, which turn into each other. Thus the 
least possible energy for this transition determines 
the ideal kinetic energy for the string dynamics. 

Since the average frequency (per step) of both ele- 
ments of a block moving upward is and both down- 
ward is (l—p)^, the average fraction of the blocks that 
change per step is -I- (1 — p)"^. This fraction times 
N/2, the total number of blocks in the string, gives 
(using Equation 17 1 the least possible average kinetic 
energy (in energy units of /i/2t). Since the motion 
exactly repeats after a cycle of N/2 steps, the ideal 
kinetic energy per cycle is invariant. 

The minimum value oi p^ + {1— p)'^ is 1/2 and oc- 
curs when p — 1/2, which is also when v = 0. Thus at 
low speeds, half of the block changes don't contribute 
to the net vertical motion of the string. If we think 
of the string as a particle moving in one dimension 
(i.e., vertically) with a horizontal internal dimension 
of width N, then only the portion of the ideal kinetic 
energy that contributes to vertical motion is particle 
kinetic energy. Thus at low speeds, the ideal parti- 
cle kinetic energy is the excess over 1 /2 of the blocks 
that change. 



N f ^ , 1 



\ + v 



1 - V 



(20) 



Since E — 7m, for small v we have E — m and we 
recover the expected non-relativistic kinetic energy. 

As p approaches zero or one, the fraction of the 
blocks that change approaches one and the vertical 
speed of the string also approaches one. Thus as w ^ 
1, all of the ideal kinetic energy contributes to the 
string motion: T ^ E. 



If we think of the block changes that don't contribute 
to overall string motion as internal dynamics of the 
string, then as the string approaches the speed 1 the 
internal dynamics stops: all of the changes contribute 
to overall string motion and none to internal dynam- 
ics. The internal dynamics exhibits relativistic time 
dilation. There is a slight subtlety, though, that 
arises from relativistically interpreting all gradient 
segments as moving. We correspondingly interpret 
a fraction 1 /7 of all segment motion as contributing 
to the internal motion, and so the internal dynamics 
slows down by this factor as the string speeds up. The 
remaining fraction 1 — 1/7 of the segment motion con- 
tributes to the overall string motion. Of course some 
of the internal "motion" we're counting here doesn't 
involve block changes, but all of the segment motion 
that contribute to overall motion of the string does 
actually correspond to block changes. Thus the frac- 
tion of all N/2 blocks of the string that contribute 
to kinetic energy of string-particle motion is in fact 
1 — 1/7 (i.e., T = E — m). The fraction of blocks 
that contribute only to internal kinetic energy in the 
string is then the total fraction that change minus 
the external kinetic fraction, and so 



N 

T]„tornal = [ + 0- ^ PY 



1-1 

7 



El'--'-^ 
.7 2 



(21) 



which is always non-negative and which approaches 
a fraction of all blocks of 1/2 — 0{v'^) as v approaches 



zero 



Related models of diffusive behavior which make 
contact with relativity are discussed in [531 13 HO] • 
None of these relativistic discussions define a model 
with objects that have an internal dynamics. 



^''Notice that in this model, for a fixed string length TV the 
energy of the string is fixed and the rest mass -E/7 approaches 
zero for strings that move at speeds approaching 1. 
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Figure 7: A collision in the Soft Sphere Model finitary 
CM dynamics. 

5 A colliding ball model 

In this section we construct a new invertible parti- 
tioning dynamics that has a large-scale limit with 
macroscopic objects and forces. This finite-state dy- 
namics is then analyzed as the integer-time behavior 
of a continuous relativistic CM dynamics. We use the 
state-change energy bound of Equation [17] to define 
the ideal energy of the model and we compare ideal 
kinetic energy with relativistic kinetic energy. 

5.1 Soft Sphere Model 

In Figure [7] we illustrate the soft sphere model 
(SSM) 0\. The SSM is an invertible and energy 
conserving CM model of classical computation simi- 
lar to Fredkin's billiard ball model. Unlike the BBM, 
however, it can be turned into a lattice gas dynam- 
ics [321 1^ with point interactions. 

As in the BBM we restrict the initial positions and 
velocities of balls to discrete values, to guarantee that 
the SSM acts like a digital system. In the collision 
shown in Figure [7] the horizontal component of veloc- 
ity of all balls entering from the left is one column per 
time step, and so consecutive moments of the history 
of the collision occur in consecutive columns. 

The collision shown is energy and momentum con- 
serving and the compression and rebound takes ex- 
actly the time needed to displace the colliding balls 
from their original paths onto the paths labeled AB. 
If a ball had come in only at A with no ball at B, it 
would leave along the path labeled AB. This model 



(a) 



(b) 



Figure 8: Rescalable soft sphere partitioning dynam- 
ics, (a) A time history of one column of particles 
colliding. At the times shown all particles are con- 
verging towards the center of the square they are in. 
(b) Site update rule. After the update particles are 
diverging from the center of the square. Only these 
cases interact (and rotations); otherwise particles go 
straight. 



is equivalent to a lattice gas automaton, with lattice 
sites located at the corners of the grid shown in the 
figure. 

5.2 Rescalable SSM 

Figure |8] describes the Rescalable Soft Sphere (RSS) 
model which implements a partitioning version of the 
SSM using colliding blocks of I's of any size. This 
is a new invertible model that conserves energy and 
momentum. Colliding blocks undergo collisions and 
then reform in appropriately shifted positions to ex- 
actly implement the SSM dynamics on any scale. The 
basic idea is that the size of the incoming blocks con- 
trols how long the interaction takes, because particles 
behave differently when they are surrounded by other 
particles than when they are alone. This allows the 
shift caused by a collision to depend on the size of 
the colliding blocks. 

Figure [8^ shows a simple collision. Two groups of 
three particles approach each other (first column on 
the left), with the top group moving right and down, 
the bottom group right and up. As in Figure [T] all 
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particles have a constant horizontal velocity and so 
consecutive discrete moments in the history of the 
collision occur in consecutive columns. 

Each square in Figure [8^ is a block of the even-step 
partition and columns are shown at consecutive even 
times. This means that all particles within a square 
are converging towards the center of the square at the 
times shown. For example, since the top three black 
particles in the first column are in upper left squares, 
they are headed down and right. To make the dia- 
gram easier to understand we've spread out the par- 
ticles in the initial state so that particles only ever 
interact on even steps, at the moments shown. The 
odd-step interaction is turned off for now — particles 
just go straight on odd steps. Later we'll make both 
steps the same and put the particles closer together. 

The size of the incoming groups of particles con- 
trols how long the interaction takes. Pairs of corre- 
sponding incoming black particles collide along the 
axis of symmetry between the incoming groups, each 
pair colliding at a separate spot. Each pair of collid- 
ing black particles turns into a pair of white particles 
moving in the same directions. The white particles 
move straight as long as they encounter other parti- 
cles, but turn back towards the axis of the collision as 
soon as they find themselves alone. The white parti- 
cles are labeled with a chirality (+ or — ) at the time 
they are created, so that they will know which way 
to turn. The white particles ignore other white par- 
ticles of the same chirality and don't interact again 
until the -I- and — come back together at the axis of 
the collision. They then turn back into black parti- 
cles and reconstitute two groups of black particles, 
each of which moves off as a unit. 

5.3 Transition rule 

The transition rule used in Figure |8^ is shown in Fig- 
ure [8|d. On the left we see the particles before the 
interaction, as they converge towards the center of 
the square. On the right we show the particles af- 
ter the interaction, as they diverge. In addition to 
the cases shown, each orthogonal rotation of a state 
on the left turns into the same rotation of the corre- 
sponding result case on the right. This allows us to 
have the same collision occur in any orthogonal oricn- 




Figure 9: A larger collision in RSS dynamics. Two 
columns of 8 particles collide. Partitions are centered 
at middles and corners of the squares shown. 

tation. In all cases not shown, there is no interaction. 
Particles go straight through the center of the square 
and continue moving in the same direction without 
change. 

This rule can be inferred from Figure [8|l. All of 
the cases shown are needed in order to reproduce 
the given evolution except for the second case (white 
to white collision) which doesn't occur. This case is 
added to make the overall dynamics invertible. Since 
cases that aren't shown don't interact, particles will 
often pass through each other without affecting each 
other. For example, if two groups of black particles 
collide head on, they pass through each other. Since 
most of the "no interaction" cases don't occur in the 
kind of collision shown in Figure [8^, we could aug- 
ment the rule without affecting Figure |8^ by changing 
some of the "no interaction" cases. 

5.4 Square balls 

Figure [9] shows a continuous time history for a col- 
lision of two groups of 8 incoming particles (black 
arrows on the left). The solid lines trace the paths 
of black particles, the dotted lines show the paths of 
white particles. Interactions occur at both the cen- 
ters and the corners of the squares. The grid is shown 
with 8x8 blocks of cells outlined by darker lines. No- 
tice that the upper group of 8 incoming black parti- 
cles occupies the fifth column within an 8x8 block, 
and the corresponding outgoing group of black par- 
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tides also occupies the fifth column within an 8x8 
block. Because of the uniform horizontal motion in 
this diagram, if we filled all of the columns in the 
two 8x8 blocks on the left of the diagram with black 
particles (moving right and down in the top block, 
right and up in the bottom block), each input col- 
umn would separately turn into the corresponding 
output column within the rightmost 8 columns. Thus 
two 8x8 "square balls" would collide, reconstituting 
themselves one ball-width to the right of where the 
balls would have gone if no collision had occurred. 
Thus the RSS dynamics implements the Soft Sphere 
collision using square balls. By making balls square 
we allow them to sometimes collide along a horizontal 
axis and sometimes along a vertical axis — this couples 
the two dimensions. 

5.5 Rescaling 

If the square balls were larger, they would shift corre- 
spondingly further to the right: the block dynamics 
is scale invariant. The square ball dynamics scales 
all the way down to 1x1 square balls: even single 
particles reproduce SSM collisions. We can take a 
Soft Sphere Model dynamics that is simulated at the 
finest grain size of the dynamics and reproduce it in 
a system m times larger in each dimension at a rate 
m times slower by simply replacing each block of the 
partition in the fine grain initial state by a tessela- 
tion of m X m identical blocks in the scaled system 
(cf. rescaling discussion in ^5]). 

Rescaling can be used as a way to take the con- 
tinuum limit of an RSS model. If the time scale of 
the underlying RSS dynamics is m times smaller than 
that of a scaled system, for m sufficiently large the 
microscopic space and time scales can be considered 
infinitesimal compared to the "macroscopic" scales. 
At that point we have the macroscopic system per- 
forming a dynamics that is exactly equivalent to a 
microscopic dynamics, but in an effectively continu- 
ous space and time. 

When macroscopic blocks collide, there is a ques- 
tion of which block is which afterwards. This of 
course depends on how we decide to interpret which 
particle is which when two particles collide. For ex- 
ample, if we decide that particle labels never cross 



the plane of a two-particle collision, then block col- 
lisions look like Figure |7] (repulsive collision) . If we 
decide particle labels always cross, then entire blocks 
pass through each other as they come together and 
then again as they come apart (attractive collision). 

5.6 ID and 3D 

In Figure [9] we analyzed the collision of 2D blocks by 
observing that each column collides independently, as 
a ID system. We could of course simply reinterpret 
the horizontal axis in the diagram as time, and this 
becomes a diagram of a ID collision of two extended 
objects. 

In this case blocks in the transition rule of Fig- 
ure should be half as wide: particles have only 
two directions (up and down). This is the only "ro- 
tation" of the rule in ID. We still need two white 
particle types so that the white particles can pass 
through other white particles of the same type un- 
til the two different types come back together and 
collide at the "axis" of the collision. 

To get the same rescalable SSM behavior in a 3D 
system we can simply apply the planar rule of Fig- 
ure [8]d in various orientations: only pairs of particles 
collide and the entire collision process takes place in 
the plane defined by the two particle directions. We 
need to use more kinds of labels for the white parti- 
cles (more possibilities of which way they can turn), 
but for each collision plane only two kinds of turn- 
ers are needed — these become the + and — particles 
in the rule. Of course the scalable 2D square balls 
become 3D cubical balls, and we only use a set of 
directions for which collisions occur between parallel 
faces of cubes. 

5.7 Mass energy 

Since the RSS model is a single-speed dynamics that 
conserves energy and momentum in its collisions, we 
are free to interpret it as a special case of a continuous 
classical relativistic dynamics with constrained initial 
conditions. 

As two particles cross over the same spot in Fig- 
ure |9] we can consider the continuation of each parti- 
cle to be either the one that went straight or the one 
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Figure 10: Macroscopic limit of RSS collision (relativistic) . (a) Fraction of energy that is mass, as a function 
of time. Three different conventions defining mass energy correspond to three different time-histories (solid- 
line, fine-dashed, coarse-dashed), (b) Vertical position y of centers of the top and bottom balls as a function 
of time, (c) Fraction of energy that is mass, as a function of separation of ball centers. 



that turned. In the former case the particles move at 
the speed of light and are massless. In the latter case 
the net motion of each particle during the timestep 
is horizontal at less than the speed of light and so it 
has mass. 

In this model all isolated freely moving particles 
have the same energy e, whose minimum possible 
value (Equation 17 1 sets the ideal energy scale. For 



a classical particle moving at the speed of light, 
p — e/c. In this section we use units in which c~\, 
so that the magnitude of the momentum equals the 
energy for a massless particle. For each particle that 
changes direction in a crossover the average momen- 
tum during the timestep is only e / ^/2 and so the mass 
is also e/\/2, since vn? — —p^. 

Thus during the collision of Figure [9] the mass goes 
from zero to a maximum of 8-\/2e and then back to 
zero. In units of v^e, the mass at each step of the 
collision is just the number of mass-generating cross- 
ings in a column. A collision of two nxn square balls 
corresponds to n simultaneous collisions such as the 
one in the figure. As we let the balls get larger and 
larger, the fraction of the total energy that is mass 
energy approaches one of the three curves shown in 
Figure [lOji. 

Here we've chosen our unit of time to be the time it 
takes a freely moving square ball to move the length 
of its diagonal. With this unit, a collision takes two 
units of time. Because of the symmetry of the col- 



lision, the fraction of the energy that is mass is the 
same for each of the two colliding balls: each ball 
changes energy into mass and then back again ac- 
cording to the curves shown. 

Which curve the mass of each ball follows depends 
on how we interpret the particle crossings. We've 
shown three possible interpretations. The "maximum 
mass" interpretation, (coarse-dashed line) is given as 
a function of time by considering all crossings to be 
massive. The "sudden mass" interpretation (thick- 
solid line) interprets all particles of each square ball 
as moving straight (with perhaps a change in particle- 
type) until the moment the two balls completely over- 
lap, and then all crossings have mass until the two 
balls start to separate. The "white mass" interpre- 
tation (fine-dashed line) is given by only counting 
pairs of white particles that cross to be massive. The 
"maximum mass" interpretation corresponds to the 
minimum possible motion (least kinetic energy). 

Let us associate the entire mass of each ball with 
a single representative point. For a freely moving 
ball, the representative point is the center of the ball. 
Since E = 7M the ratio of mass to energy shown in 



Figure 10 J, is equal to I/7. Knowing 7 as a function 
of time gives us the velocity v of the ball as a func- 
tion of time. Since the horizontal component of v is 
unaffected by the collision and remains constant this 
gives us the vertical velocity as a function of time, 
and hence the vertical position as a function of time. 
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This is depicted in Figure 10 d, where we've shown 
the path that the representative point for each ball 
must follow during a collision if its mass varies ac- 
cording to one of the three curves in Figure [TO^. The 
three cases show that exactly the same collision can 
be interpreted as attractive, repulsive, or sticky. 

It is clear from the behavior of M/E and y as 
functions of time that the ratio M/E depends only 
on the separation r between the two representative 
points (once the interaction starts) . This dependence 
is graphed in Figure 10". Thus if we think of mass 
as a kind of relativistic potential energy (i.e., total 
energy minus kinetic energy), we see that it depends 
only on the distance between the two balls. The po- 
tential is zero when the balls are not yet close enough 
to touch and it increases as the balls collide. The 
maximum mass M for each ball is the fraction E j\f2 
of the ball's energy, which obtains when all particles 
are paired so that there is maximum cancellation of 
vertical momentum components. 



tion 17 1. In constrast, if the relativistic kinetic en- 



ergy of a free massless particle is e, then the relativis- 
tic kinetic energy of a block containing two massless 
particles whose paths intersect at right angles is only 
2e(l - 1/7) « .6e. 

We can try to reinterpret the RSS dynamics to 
make the relativistic kinetic energy match the state- 
change energy both for a free particle and for a two- 
particle collision. If we let free RSS particles each 
have mass m/ and let a two-particle collision have 
mass Trie we find that the only choice of masses that 
allows the two energies to be equal is mc = 2m j, 
which corresponds to 7 = 1. Thus the state-change 
energy becomes the ideal kinetic energy in the RSS 
model only for a non-relativistic interpretation of the 
dynamics. 



5.8 Relativistic invariance? 

We were able to use the RSS model to discuss mass in 
a relativistic collision by interpreting the model's sin- 
gle particle-speed as the speed of light and its conser- 
vations as relativistic conservations. This model is, 
however, clearly not very relativistic even in a con- 
tinuum limit. While it does display macroscopic me- 
chanical behavior, it does so only in a single inertial 
frame. It thus does not exhibit relativistic frame in- 
variance. 

We might expect that, for a colliding pair of RSS 
particles, the relativistic kinetic energy (which in- 
cludes only energy related to net translational mo- 
tion) should be closely related to the state-change 
energy (which is the minimum energy for the net mo- 
tion) . This suggests that we might be able to use the 
intrinsic definition of ideal kinetic energy provided by 
state-change to expose the non-relativistically invari- 
ant character of the RSS model. 

If wc interpret freely moving RSS particles as mass- 
less, as we did in Section |5.7| then state-change en- 
ergy and relativistic kinetic energy evolve differently. 
The state-change energy of both a free particle and of 
a block where two particles cross are the same (Equa- 



6 Conclusions 

A finite-sized QM system has only a finite number 
of possible distinct states and can move between dis- 
tinct states at only a finite rate. This means that 
classical special cases of QM dynamics must be clas- 
sical finite-state dynamics. Conversely, we can regard 
the dynamics of real physical systems to be a general- 
ization of classical finite-state dynamics. This is not 
a common viewpoint among physicists today. 

The study of classical finite-state dynamics that 
are special cases of QM dynamics should, at the least, 
be of interest for pedagogical reasons since it allows 
physical concepts to be seen in an intuitive classi- 
cal setting. This study should also be of interest 
to theoretical computer science since it establishes 
an exact correspondence between a classical compu- 
tation and an ideal quantum realization. Finally, 
this study should be of interest to theoretical physics 
since it provides a novel finite-state perspective on 
the foundations of both classical and quantum me- 
chanics, along with an intuitive starting point for the 
construction of new physical models. 
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